Polymers in a vacuum 
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In a variety of situations, isolated polymer molecules are found in a vacuum and here we examine 
their properties. Angular momentum conservation is shown to significantly alter the average size of 
a chain and its conservation is only broken slowly by thermal radiation. The time autocorrelation for 
monomer position oscillates with a characteristic time proportional to chain length. The oscillations 
and damping are analyzed in detail. Short range repulsive interactions suppress oscillations and 
speed up relaxation but stretched chains still show damped oscillatory time correlations. 



The properties of polymer chains have been investi- 
gated extensively over the past fifty years [1] but the 
vast majority of these studies have been concerned with 
situations where they are in a solution or a melt. How- 
ever there are some situations where polymer molecules 
are essentially in a vacuum. Desorption and ionization 
of polymers, often by lasers is carried out during mass 
spectrometry, in order to characterize desorbed proteins. 
This has many important applications including the un- 
derstanding of cancer [2] . Polymer molecules of many dif- 
ferent kinds have been detected in interstellar media [3] 
and although they have predominantly had less than 10 
units, the detection of new species is an active area of 
research. It might also be possible to employ optical 
tweezers on biomolecules such as DNA and manipulate 
them in a vacuum, in a manner similar to what is now 
done routinely in aqueous solution [4]. 

To aid in the possible experimental observation of such 
systems, some basic properties of isolated polymers in 
a vacuum are considered here. The first question that 
we ask is how their statistics are modified from those 
in solution. Solvents will compete with intra-chain at- 
tractions so that above the temperature [1], a polymer 
chain will be swollen. Without the solvent present, this 
would imply that a chain at the same temperature would 
be collapsed. But at high enough temperatures, entropy 
will dominate over energy, and a polymer, just like a liq- 
uid, will then want to expand into a gas, or self avoiding 
phases. Because carbon-carbon bonds are very strong, 
it might then be possible to find some species where a 
polymer will become swollen in isolation for long enough 
periods of time to be observable. Even if it turns out that 
this is not possible, polymers through desorption often 
carry a charge, and this additional coulomb repulsion is 
quite substantial, at 500K it is ~ SSfc^T for two electrons 
1 nm apart. This will serve to stretch a chain. Such a 
situation is a possibility in the desorption and ionization 
of proteins done in mass spectrometry experiments [2]. 

It might then appear that the statistics of this sys- 
tem are identical to that of a chain in a solvent, with 
some modification of interaction parameters. However 
one important difference is the conservation of angular 
momentum that we might expect to see in this case as 
opposed to a polymer in a solvent. In statistical mechan- 



ics, this conservation law is ordinarily ignored and is not 
expected to make a difference to system properties when 
the number degrees of freedom are large. However we 
will sec that for a polymer in isolation it has a significant 
effect on its size, even when the total angular momentum 
is zero. The effects of angular momentum conservation 
have been recently studied in self-gravitating systems [5] 
where it leads to different phases for some models for 
finite angular momentum. 

The starting point for this situation is the formula for 
the classical entropy in the microcanonical ensemble with 
conservation of linear and angular momentum of a system 
with total potential energy U [5]. However we can safely 
transform this into the canonical ensemble at tempera- 
ture T using the usual argument that the fluctuations 
in the energy at constant temperature are small for a 
large number of degrees of freedom. However we are still 
keeping the S function constraint on the total angular 
momentum L. Using the Fourier representation for a S 
function and integrating over the momenta, the partition 
function becomes 
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where N is the number of monomers, I the moment of 
inertia tensor of a polymer conformation, and the center 
of mass is fixed to zero. Adding a term e J^i '''i to /3U 
allows one to differentiate In Z with respect to e to obtain 
the r.m.s. size of a chain defined as (i?^) = -^(X^i^i 

For an "ideal" or "phantom" chain [6] with ring topol- 
ogy, the calculation can be done exactly [7] giving the 
results shown in Fig. 1. The rescaled angular momen- 
tum L' = LVl2/{NVTml) where m is the mass and I is 
the step length, which both can be taken equal to 1. 

At L = 0, {R^)/N — .07. Results using a simulation 
method described below give .071 which are the same 
to within the error bars. This is substantially below, 
1/12 « .083, the value of the same quantity where con- 
servation of angular momentum is not enforced. When 
averaged over all angular momenta, the size of the chain 
must agree with the non-conserved case, and because 
high L chains will have greatly extended conformations. 
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FIG. 1; (Color Online) The variation of the size of a chain 
versus its angular momentum for a ideal ring chain. The 
straight line is the asymptotic slope which can be obtained 
through a simple analysis. 



this must be compensated for by correspondingly com- 
pact configurations for small L. 

To obtain the asymptotic behavior for large L, using 
a much simpler argument, we expect that in this limit 
the dominant configuration of the chain will be a highly 
stretched circle of radius i? rotating symmetrically about 
the axis of angular momentum. We minimize the free en- 
ergy, of a polymer taking into account both its kinetic en- 
ergy and elastic energy yielding /N = L'/(127r) again 
choosing I ~ m ^ 1. The straight line in fig. 1 has the 
same slope. 

The case of self avoiding (swollen) chains {R ~ N'' 
with u « 3/5), is qualitatively similar, and the asymp- 
totic scaling can be easily worked out along similar lines 
as above, giving R oc L^^^^N^^^^ . 

The total angular momentum however, is not con- 
served. Interaction with thermal photons will cause the 
angular momentum to equilibrate on a time scale that 
we will now estimate. First we consider the flux of elec- 
tromagnetic energy emitted by a single polymer. The 
emission of thermal radiation per unit area of a black 
body is given by the Stephan-Boltzmann law S — aT'^, 
where a is the Stephan-Boltzmann constant. However 
this greatly overestimates the radiation because of the 
weak efficiency of small objects in emitting light of a far 
greater wavelength. Calculations for metal nanoparticles 
(which should be better emitters than dielectrics) give a 
suppression factor of less than 10~^ when the nanopar- 
ticles are 50 A in radius [8]. This gives a ratio of ksT 
to emitted power of 6 x IO~^s. The relaxation time as 
we will see below scales roughly as N'^. The microscopic 
hopping time typical for such a system is 10~"'^^s, meaning 
that the relaxation time for N = 100 is roughly 10~^s. 
Thus thermal photon equilibration is more than two or- 
ders of magnitude slower than the time-scale for relax- 
ation of a chain. Therefore one expects to see transitions 



in the time averaged radius of gyration of a chain as 
photons are emitted and absorbed by the polymer. This 
might be observable in the signature of noise seen in light 
scattering. 

We now turn to a study of the dynamics of these poly- 
mers and we will see that in this respect the situation 
is very different from that of a polymer in a solvent. In 
a solvent, models for dynamics for the most part con- 
sider monomers connected together by linear springs, 
thermal noise, solvent drag, and hydrodynamic interac- 
tions. However in a vacuum, the last three terms are not 
present, which leaves us with Newtonian dynamics for 
linear springs, but this will never thermalize, so nonlin- 
ear forces must be considered. 

One might first guess that sufficient nonlinearity would 
introduce strong enough dissipation of individual modes 
so that the behavior would be similar to that of the 
Rouse model [9], which describes a "free draining" chain. 
That is one where individual monomers experience a 
drag proportional to their velocity. However this violates 
Galilean invariance, as we shall discuss in more detail 
below. Moreover one dimensional systems are notorious 
for not being able to equilibrate energy well, and there 
are many well studied instances where this problem is 
known to occur. Perhaps the best known example of this 
is the Fermi Pasta Ulam chain [10]. For small enough en- 
ergies this system shows strong recurrences in amplitude 
of modes above which it equilibrates [11] Models that 
are even more nonlinear such as the Sinai-Chernov "Pen 
Case model" [12] do not suffer from this problem, how- 
ever they still exhibit highly non-local time correlations, 
with universal power law decays [13], which is a general 
result for one dimensional chains that are momentum and 
energy conserving [14] 

So we first consider the dynamics of a chain neglect- 
ing any self avoiding interactions, but using an athermal 
highly nonlinear model for the reasons mentioned above. 
Thus we have chosen a model where monomers of equal 
mass are coupled together by links of fixed length. Aside 
from this constraint, there is no potential energy. The 
monomers can freely rotate but there is no coupling to an 
outside system so that there is no dissipation or random 
noise term. The model rigorously satisfies conservation 
of energy, momentum and angular momentum. An effi- 
cient method for evolving such chains was developed so 
that despite the large number of length constraints, the 
computation for each time step scales linearly with the 
number of monomers. The details will be published else- 
where [7]. The angular momentum, center of mass and 
energy were monitored to ensure that their drifts due to 
numerical error remained small for all data used. 

The monomer-monomer autocorrelation function de- 
fined as 
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FIG. 2: (Color Online) The autocorrelation function g{t) as 
defined in Eq. (2) for a chain of 128 monomers. The inset 
shows the beginning of the plot at a higher magnification, 
where the initial linear increase and cusps are apparent. 



was calculated for chains of different lengths and is dis- 
played in Fig. 2 for N = 128 averaged over 24, 000 runs. 
This is very unlike the correlation function for a polymer 
in solution which shows a smooth slow increase, not the 
wildly oscillatory form seen here. The period of these os- 
cillations scales as chain length. This is what one would 
expect for a linear model with no friction, because the 
lowest mode of oscillation has spring constant oc 1/A^, 
and the mass is cx iV giving a period cx N. At first 
sight, it is surprising that the oscillations are so weakly 
damped in a model that is so highly nonlinear. To un- 
derstand this better, we analyze this problem in terms 
of linear modes. The relaxation time of the chain is 
then given by the damping time for these oscillations. 
In terms of a Fourier decomposition along the arclength 
of the chain, we can model the correlation function as 
(rfc(G)r^(i)) ^ {\rk\'^) Re{cxp{{iujk - Afe)<) and because 
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and this can be used to fit the numerical data. In this 
expression the w^'s represent the frequencies of oscilla- 
tion. We therefore expect that for small k that Afc will 
be small in comparison to ujk, the only other relevant 
frequency scale for each term in the summation. 

For small t, and no dissipation (Afc = 0), the small 
t behavior of g{t) in the above expression can easily be 
shown to be oc t. Yet for short times this shows that 
the mean square displacement follows the same law as a 
diffusive process. For longer times, oscillatory behavior 
will be seen but will be asymmetric, showing cusps at the 
minima and parabolic maximum. These same qualitative 
features persist for small and are also in quite good 



agreement with the simulation data shown in Fig. 2. 

It is clear from the data that LUk oc fc, just as one would 
expect from a linear system of springs. The damping ap- 
pears to fit best to a form close to Afc oc A:^^ Fitting this 
to different chain lengths, N = 64 and 128, gives a relax- 
ation time Trei oc ^('^■^^±.15) _ ]^qIq n^g^^ ^j^g ^ase of 
one dimensional heat conduction, it has been found that 
even with highly nonlinear models [13], asymptotic large 
N behavior is difficult to study as more than 10^ particles 
must be considered to get a good estimate of critical ex- 
ponents. Therefore it is possible that the exponent found 
is off by ^ 10% of its asymptotic value. 

This value of the relaxation time scales closely to what 
is seen in the Rouse model although in that case the 
physics is very different as there is no inertia term and 
no rapidly oscillatory behavior. In the case we consider 
here, the origin of this time is due to nonlinear coupling 
of different modes resulting in the slow translation to 
decoherence as time progresses. 

As mentioned above, this problem is quite similar to 
that of a one dimensional nonlinear chain of particles. In 
that case the system is also characterized by long wave- 
length oscillations that slowly decohere. But there, the 
relaxation time Trei oc N^/^ which is different than the 
polymer case. If the polymer chain was stretched by a 
constant force so that it was quasi-one dimensional, one 
would expect the same iV'^/^ scaling for the relaxation 
time. 

To understand the damped oscillatory behavior of the 
correlation function, in Eq. 3, we consider what kind 
of coarse grained linear stochastic equation would best 
approximate the evolution of r(s,i), the position of the 
chain at arclength s and time t. Because the system has 
Galilean invariance, there can be no r term for the damp- 
ing, as center of mass velocity is conserved. By symme- 
try, the lowest order damping term must be Cd^r/dtd^s 
where C is a constant. Adding in inertia, random forcing 
^(s,t) and chain connectivity gives 
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In relation to the above analysis, this gives a damping 
Afc cx for small k. Such a model matches fairly well 
the numerical data [15]. 

Such a proposal for internal damping for polymer 
chains in connection with Cerf friction [1, 16] has been 
made before [17] using a nonrigorous derivation that gives 
rise to the same third order derivative term as in eqn. 4. 
Adding such a term to the Rouse equation provides an 
explanation of experiments [18] on extensional relaxation 
of polymers in solvent. Solvents with different viscosities 
were considered to extrapolate to the limit of zero vis- 
cosity, and the results can be interpreted using such a 
term. 

As one might expect, the inclusion of repulsive interac- 
tions between monomers suppresses the oscillations that 
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FIG. 3: (Color Online) The scaled autocorrelation function 
g{t)l{R^) versus t/W, for three chain lengths, = 32, 64, 
and 128 with short-range repulsive interactions. The inset 
shows g{t) for an A'^ = 32 chain with equal charges at both 
ends that cause it to stretch. 



rameters were chosen so that there is still a substantial 
amount of interaction between neighboring monomers, 
but the chain is still quite stretched. Here one can clearly 
sec oscillations in g{t)^ intermediate in behavior between 
the ideal chain and interacting cases. 

In conclusion, the equilibrium statistics and dynamics 
of polymers in a vacuum have many interesting prop- 
erties. The addition of angular momentum conservation 
significantly alters chain statistics. The subtle power-law 
time correlations found in momentum conserving one di- 
mensional systems can lead to dynamics that are oscil- 
latory and show unusual scaling properties. It is hoped 
that this work will provide impetus for further experi- 
mental observation of these fascinating systems. 

The author thanks Larry Sorensen for very useful dis- 
cussions. 



are seen in the ideal chain. For the sake of efficiency, soft- 
core potentials were added between between monomers 
having a potential of the form l3V{r) = 2(1 — (r//)^)^. 
Statistics of such chains with total angular momentum 
of zero were measured and the size scaling exponent gave 
v — .596 ± .01 in good agreement with the well known 
three dimensional value. 

In Fig. 3 the autocorrelation function for these chains, 
see Eq. 2, is plotted for chain lengths N = 32, 64, 
and 128, on scaled axes so that they coincide for short 
times. The vertical axis is g{t)/{B?) and the horizon- 
tal one is t/NP, with p chosen to fit short times best. 
With p = 1.15, the plots coincide well over half of the 
vertical range, from to 1. However the long time be- 
havior for iV = 32 is noticeably above the longer length 
chains. However = 64 is only slightly above N — 128, 
and given the correlated error bars, this is barely statis- 
tically significant. This is strong evidence that for large 
N the correlation function approaches the scaling form 
g{t) — N'^'^ f {t / NP) , and therefore the relaxation time for 
this chain is oc Np, with p = 1.15 ± .05. Note that this 
is much smaller than that of the ideal chain discussed 
above, presumably because long range interactions along 
the chain backbone allow much faster equilibration of en- 
ergy and momentum [19]. We expect the time it takes a 
chain segment to move of order its average size Rg should 
be Rg divided by the center of mass speed for of order 
half the chain, ~ N~^/^, which gives trei ~ N^-^. 

Finally we contrast this with what happens if charges 
are added to both ends. With charged protein molecules 
observed in mass spectrometry, a similar situation could 
also occur. The inset in Fig. 3 shows the autocorrelation 
function for this case, where the end to end distance is 
10.03, about one third of the chain's arclength. The pa- 
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